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There is growing interest to investigate states of matter with topological order, which support excitations in the 
form of anyons, and which underly topological quantum computing. Examples of such systems include lattice 
spin models in two dimensions. Here we show that relevant Hamiltonians can be systematically engineered 
with polar molècules stored in optical lattices, where the spin is represented by a single electron outside a closed 
shell of a heteronuclear molecule in its rotational ground state. Combining microwave excitation with the 
dipole-dipole interactions and spin-rotation couplings allows us to build a complete toolbox for effective two- 
spin interactions with designable range and spatial anisotropy, and with coupling strengths significantly larger 
than relevant decoherence rates. As an illustration we discuss two models: a 2D square lattice with an energy 
gap providing for protected quantum memory, and another on stacked triangular lattices leading to topological 
quantum computing. 



I. INTRODUCTION 

Lattice spin models are ubiquitous in condensed matter 
physics where they are used as simplified models to describe 
the characteristic behavior of more complicated interacting 
physical systems. Recently there have been exciting theoret- 
ical discoveries of models with quasi-local spin interactions 
with emergent topological order [1,2]. In contrast to Landau 
theory where various phases of matter are described by bro- 
ken symmetries, topological ordered states are distinguished 
by homology class and have the property of being robust to 
arbitrary perturbations of the underlying Hamiltonian. These 
states do not exhibit long range order in pairwise operators, 
rather they have long range order in highly nonlocal strings 
of operators. A real world example is the fractional quantum 
Hall effect which gives rise to states with the same symmetry 
but distinguishable by quantum numbers associated with the 
topology of the surface they live on [3]. 

It is of significant interest to "design" materials with these 
properties, both to observe and study exòtic phases, and in 
light of possible applications. Cold atòmic and molecular 
gases in optical lattices are prime candidates for this endeavor 
in view of the complete controllability of these systems in 
the laboratory. The idea of realizing bosonic and fermionic 
Hubbard models, and thus also lattice spin models, with cold 
àtoms in optical lattices has sparked a remarkable series of 
experiments, and has triggered numerous theoretical studies 
to develop cold àtoms as a quantum simulator for strongly 
correlated condensed matter systems (see e.g. [4, 5, 6] and 
references therein). However, coaxing a physical system to 
mimic the required interactions for relevant lattice spin mod- 
els, which mustbe both anisotropic in space and in the spin de- 
grees of freedom, and given range, is highly nontrivial. Here 
we show that cold gases of polar molècules, as presently de- 
veloped in the laboratory [ ], allow us to construct in a natural 
way a complete toolbox for any permutation symmetric two 
spin- 1/2 (qubit) interaction. The attractiveness of this idea 
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also rests on the fact that dipolar interactions have coupling 
strengths significantly larger than those of the atòmic Hubbard 
models, and relevant decoherence rates. 

Our bàsic building block is a system of two polar molècules 
strongly trapped a given sites of an optical lattice, where the 
spin- 1/2 (or qubit) is represented by a single electron out- 
side a closed shell of a heteronuclear molecule in its rotational 
ground state. Heteronuclear molècules have large permanent 
elèctric dipole moments. This implies that the rotational mo- 
tion of molècules is coupled strongly via the dipole-dipole in- 
teractions, whose signatures are the long range 1 /r 3 character 
and an angular dependence, where the polar molècules attract 
or repèl each other depending on the relative orientation of 
their dipole moments. In addition, microwave excitation of 
rotational energy levels allows to effectively tailor the spatial 
dependence of dipole-dipole interactions. Finally, account- 
ing for the spin-rotation splitting of molecular rotational levels 
we can make these dipole-dipole interactions spin-dependent. 
General lattice spin models are readily built from these binary 
interactions. 



II. ANISOTROPIC SPIN MODELS WITH NOISE 
RESILIENT GROUND STATES 

Two highly anisotropic models with spin- 1/2 partides 
which we will show how to simulate are illustrated in Figs. la 
and lb respectively. The first takes place on a square 2D lat- 
tice with nearest neighbor interactions 

= t't'Wj+i +cosÇo^o^ w ). (1) 

'=1 7=1 

Introduced by Duoçot et al. [8] in the context of Josephson 
junction arrays, this model (for Ç ^ ±7t/2) admits a 2- fold 
degenerate ground subspace that is immune to local noise up 
to £th order and hence is a good candidate for storing a pro- 
tected qubit. 

The second, occurs on a bipartite lattice constructed with 
two 2D triangular lattices, one shifted and stacked on top of 
the other. The interactions are indicated by nearest neighbor 
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FIG. 1 : (Color online) Example anisotropic spin models that can be 
simulated with polar molècules trapped in optical lattices. (a) Square 
lattice in 2D with nearest neighbor orientation dependent Ising inter- 
actions along x and z- Effective interactions between the spins Si and 
S2 of the molècules in their rovibrational ground states are generated 
with a microwave field E(f) inducing dipole-dipole interactions be- 
tween the molècules with dipole moments Dj and D2, respectively. 
(b) Two staggered triangular lattices with nearest neighbors oriented 
along orthogonal triads. The interactions depend on the orientation 
of the links with respect to the elèctric field. (Dashed lines are in- 
cluded for perspective.) 



links along the x,y and z directions in real space: 
H™=J X £ <? j( ? k +J x £ a]a) • ./. £ a)a\. (2) 
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This model has the same spin dependence and nearest neigh- 
bor graph as the model on a honeycomb lattice introduced by 
Kitaev [9]. He has shown that by tuning the ratio of interac- 
tion strengths \ J± |/| J z \ one can tune the system from a gapped 
phase carrying Abelian anyonic excitations to a gapless phase 
which in the presence of a magnètic field becomes gapped 
with non- Abelian excitations. In theregime ^ 1 the 

Hamilonian can be mapped to a model with four body oper- 
ators on a square lattice with ground states that encode topo- 
logically protected quantum memory [10]. One proposal [11] 
describes how to use trapped àtoms in spin dependent optical 

lattices to simulate the spin model H^ n . There the induced 
spin couplings are obtained via spin dependent col·lisions in 
second order tunneling processes. Larger coupling strengths 
are desirable. In both spin models (I and II) above, the signs 
of the interactions are irrelevant although we will be able to 
tune the signs if needed. 



III. SPECTROSCOPY OF POLAR MOLÈCULES IN 
OPTICAL LATTICES 

Our system is comprised of heteronuclear molècules with 
2 £i/2 ground electrònic states, corresponding for example to 
alkaline-earth monohalogenides with a single electron outside 
a closed shell. We adopt a model molecule where the ro- 
tational excitations are described by the Hamiltonian H m = 
BN 2 + yN • S with N the dimensionless orbital angular mo- 
mentum of the nuclei, and S the dimensionless electrònic spin 
(assumed to be S — 1/2 in the following). Here B denotes 
the rotational constant and y is the spin-rotation coupling con- 
stant, where a typical B is a few tens of GHz, and y in the 



hundred MHz regime. The coupled basis of a single molecule 
i corresponding to the eigenbasis of H' m is {\Ni,Si,Ji;Mj t )} 
where J, = N; + S; with eigenvalues E(N = 0, 1 /2, 1 /2) 
()./:: 1.1 2 : 1 2: =2B-y, and £(1,1/2,3/2) = 2B + y/2. 
While we ignore hyperfine interactions in the present work, 
our discussion below is readily extended to include hyperfine 
effects, which offer extensions to spin systems S > 1/2. 

The Hamiltonian describing the internal and external dy- 
namics of a pair of molècules trapped in wells of an optical 
lattice is denoted by H = +H SX , The interaction describing 
the internal degrees of freedom is H m = Hdd + L^=i H' m . Here 
i/dd is the dipole-dipole interaction given below in Eq. (3). 
The Hamiltonian describing the external, or motional, degrees 
of freedom is H ex = Y%=\ P? /(2'«) + V/(xj — x,), where P, is 
the momentum of molecule ; with mass m, and the potential 
generated by the optical lattice V, (x — Xj) describes an exter- 
nal confinement of molecule i about a local minimum x, with 
1D rms width zo- We assume isotropic traps that are approx- 
imately harmònic near the trap minimum with a vibrational 
spacing /iCOosc- Furthermore, we assume that the molècules 
can be prepared in the motional ground state of each local 
potential using dissipative electromagnètic pumping [12]. It 
is convenient to define the quantization axis z along the axis 
connecting the two molècules, X2 — xi = Azz with Az corre- 
sponding to a múltiple of the lattice spacing. 

The near field dipole-dipole interaction between two 
molècules separated by r = xi — X2 is 



Hdd 



£ ((-l) q D\ D 2 - q -3D\ Q D 2ü +h.c). (3) 
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The dipole operator coupling the ground and first rotational 
states of molecule is Dt = E«=-i \ N = l>q}u{N = °,0|é 
and d is the dimensionful dipole moment. 

While the present situation of dipole-dipole coupling of ro- 
tationally excited polar molècules is reminiscent of the dipole- 
dipole interaction of between electronically excited atom pairs 
[13], there are important differences. First, unlike the atòmic 
case where electronically excited states typically are anti- 
trapped by an optical lattice, here both ground and excited 
rotational states are trapped by an essentially identical po- 
tential. Hence motional decoherence due to spin dependent 
dipole-dipole forces is strongly suppressed by the large vibra- 
tional energy Tico osc . Second, scattering rates are drastically 
reduced. The decay rate at room temperature from excited ro- 
tational states is ~ 10~ 3 Hz [14] versus a comparable rate of 
MHz for excited electrònic states. 

The ground subspace of each molecule is isomorphic to a 
spin 1 /2 particle. Our goal is to obtain an effective spin-spin 
interaction between two neighboring molècules. Static spin- 
spin interactions due to spin-rotation and dipole-dipole cou- 
plings do exist but are very small in our model: H v dw( r ) = 

-(d 4 /2Br 6 ) \l + (j/4B) 2 (1 +4Si • S 2 /3 - 25f5|)l . The first 

term is the familiar van der Waals 1 / r 6 interaction, while 
the spin dependent piece is strongly suppressed as y/4B ss 
10~ 3 <C 1. Therefore, we propose dynamical mixing with 
dipole-dipole coupled excited states using a microwave field. 
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The molècules are assumed trapped with a separation 
Az ~ ry = (2c/ 2 /y) 1 / 3 , where the dipole dipole interaction is 
^ 2 / r Y = Y/2- In this regime the rotation of the molècules is 
strongly coupled to the spin and the excited states are de- 
scribed by Hunds case (c) states in analogy to the dipole- 
dipole coupled excited electrònic states of two àtoms with 
fine-structure. The ground states are essentially spin indepen- 
dent. In the subspace of one rotational quantum (N\ +N2 = 1), 
there are 24 eigenstates of //;„ which are linear superpositions 
of two electron spin states and properly symmetrized rota- 
tional states of the two molècules. There are several sym- 
metries that reduce H m to block diagonal form. First, H^d, 
conserves the quantum number Y — Mn + M$ where Mn = 
Mjv, + M Nl _ and M s = M Sl + M$ 2 are the total rotational and 
spin projections along the intermolecular axis. Second, par- 
ity, defined as the interchange of the two molècules followed 
by parity though the center of each molecule, is conserved. 
The o = ± 1 eigenvalues of parity are conventionally denoted 
g(u) for gerade(imgerade). Finally, there is a symmetry as- 
sociated with reflection R of all electrònic and rotational co- 
ordinates through a plane containing the intermolecular axis. 
For \Y\ > all eigenstates are even under/? but for states with 
zero angular momentum projection there are ± 1 eigenstates of 
R. The 16 distinct eigenvalues correspond to degenerate sub- 
spaces labeled |T|„ (J) with J indicating the quantum number 
in the r — > °° asymptotic manifold (N = 0,/ = 1 /2;N = Í,J). 
Remarkably, the eigenvalues and eigenstates can be comput- 
ing analytically yielding the well known Movre-Pichler po- 
tentials [15] plotted in Fig. 2. 



IV. ENGINEERING SPIN-SPIN INTERACTIONS 

In order to induce strong dipole-dipole coupling we intro- 
duce a microwave field E(x,t)ep with a frequency cop tuned 
near resonance with the N = — > N = 1 transition. Be- 
cause the rotational states are spaced nonlinearly, this transi- 
tion is resolvable without coupling to higher rotational states 
by multiphoton processes. In the rotating wave approxima- 
tion, the molecule-field interaction is H m f = — £? =1 (MÏI)t • 
e F e'( kf " x ·~ (í>F ' s> /2 + h.c), where the Rabi frequency is |Í2| = 
d\Eo\/h. As the molècules are trapped by optical wavelengths 
such that kfAz ~ 10~ 5 the dipoles are excited in phase only. 

The effective Hamiltonian acting on the ground states is ob- 
tained in second order perturbation theory as 



Heff(r)=££ 

iJHr) 



(g f \H mí \X(r))(X(r)\H mí \ 8i ) 
h(0 F -E(X(r)) 



\gf)(gil (4) 



where {|gí),|g/)} are ground states with N\ = N2 = and 
{|A,(r))} are excited eigenstates of //j n with N\ +N2 = 1 and 
with excitation energies {E(k(r))}. The reduced interaction 
in the subspace of the spin degrees of freedom is then ob- 
tained by tracing over the motional degrees of freedom. For 
molècules trapped in the ground motional states of isotropic 
harmònic wells with rms width zo the wave function is separa- 
ble in center of mass and relative coordinates with the relative 




FIG. 2: (Color online) Movre-Pichler potentials for a pair of 
molècules as a function of their separation r. The potentials £(g, (r)) 
for the 4 ground-state (dashed lines) and the potentials E(k(r)) for 
the first 24 excited states (sòlid lines). The symmetries \Y\a of the 
corresponding excited manifolds are indicated, as are the asymptotic 
manifolds (Nj,Jj;Nj,Jj). The relative coordinate probability densi- 
ties on a square lattice are depicted on the ground state potential. 



coordinate wavefunction 



7t 3 /4(2 Z0 ) 3 /2 



-(r 2 +A^-2rAzcos9)/8zg 



where cos = r • z/r. The effective spin-spin Hamiltonian is 
then // spin = (Z/ e ff(r)) re i. 

The Hamiltonian in Eq. (4) is guaranteed to yield some en- 
tangling interaction for appropriate choice of field parameters 
but it is desirable to have a systematic way to design a spin- 
spin interaction. Fortunately, the model presented here pos- 
sesses sufficient structure to achieve this essentially analyti- 
cally. The effective Hamiltonian on molècules 1 and 2 induced 
by a microwave field is 



H eS {r) 



h\Q. 



3 

L < 

a,p=0 



l^a 



(5) 



where {o a } 3 x= o = {ljO^G-^G''} and A is a real symmetric 
tensor. See App. A for an explicit form of the matrix coeffi- 
cients as a function of field polarization and frequency. 

Eq. (5) describes a genèric permutation symmetric two 
qubit Hamiltonian. The components Aq s describe a pseudo 
magnètic field which acts locally on each spin and the compo- 
nents A SJ describe two qubit coupling. The pseudo magnètic 
field is zero if the microwave field is linearly polarized but 
a real magnètic field could be used to tune local interactions 
and, given a large enough gradient, could break the permuta- 
tion invariance of // sp ; n . 

For a given field polarization, tuning the frequency near an 
excited state induces a particular spin pattern on the ground 
states. These patterns change as the frequency is tuned though 
múltiple resonances at a fixed intermolecular separation. In 
Table I it is shown how to simulate the Ising and Heisen- 
berg interactions in this way. Using several fields that are 
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TABLE I: Some spin patterns that result from Eq. (5). The field 
polarization is given with respect to the intermolecular axis z and 
the frequency (üf is chosen to be near resonant with the indicated 
excited state potential at the internuclear separation Az. The sign of 
the interaction will depend on whether the frequency is tuned above 
or below resonance. 



Polarization 


Resonance 




Spin pattern 
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2 g 




a z a z 
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°J 




cr r cr v + a y a y - a z a z 
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o« 




c^'cr 1 — a y a y + a z a z 


y 


ot 




— a^a 1 + a y a y + a z a z 


{$-&)/ y/2 


°i 




—a x o y — + a z a z 


cos + sin 


h 


h 


((fcf + a z à x ) + X2<5 z a z 
+l 3 (c x c x + c y cy) 


cos ty + sin í;£ 


h 


h 


(a y a z + a'a- v ) + Xxa z a z 
+X 3 (cr t cr t + a v o v ) 



sufficiently separated in frequency, the resulting effective in- 
teractions are additive creating a spin texture on the ground 
states. The anisotropic spin model Hxyz = A, A C5 A C7 jL + X y O y O y + 
L·a z o z can be simulated using three fields: one polarized 
along z tuned to 0+(3/2), one polarized along y tuned to 
0~(3/2) and one polarized along y tuned to 0+(l/2). The 
strengths Xj can be tuned by adjusting the Rabi freqüències 
and detunings of the three fields. Using an external magnètic 
field and six microwave fields with, for example, freqüències 
and polarizations corresponding to the last six spin patterns in 
Table I, arbitrary permutation symmetric two qubit interaction 
are possible. 

The effective spin-spin interaction along a different 
intermolecular axis z! can be obtained by a frame 
transformation in the spherical basis. Writing ï = 
Z) lt (Pi,p 2 ,P3)·(0,l,0) r , where D j is the spin-j Wigner ro- 
tation, the effective Hamiltonian along £ in the original 
coordinate system is obtained by the following replace- 
ments to the field polarization vector and spin operators: 
(oL,a ,a + ) 7 ' -> D lt (Pi,p 2 ,p3).(a-,a ,a+) 7 ' and o" -> 
£) 1 / 2 (Pi,P 2 ,P3)a a £) 1 /2t(p l5 p 2 p 3 ) For example, using a z 
polarized field tuned near 0+(3/2) and a field polarized in the 
x — y plane tuned near l„(3/2) creates a Heisenberg interac- 
tion between any two molècules separated by r with arbitrary 
orientation in space. 



V. APPLICATIONS 

We now show how to engineer the spin model I. Consider a 
system of trapped molècules in a square lattice with site coor- 
dinates in the z — x plane {x/y} = {ibz + jbx; i,j G [l,í]n^}- 
Illuminate the system with a microwave field with linear po- 
larization e/r = cosÇy + sinÇi and field frequency cop tuned 
such that the peak of the relative coordinate wavefunction at 
r = b is near resonant with the 2 g potential but far detuned 



from other excited states. Then the dominate interaction be- 
tween nearest neighbor molècules is of Ising type along each 

axis and we realize with J = (h\ü.\) 2 (l /8{h& F -2B- 

y/2 — í/ 2 /r 3 )) re i. For realistic parameters, this coupling can 
range from 10 — 100 kHz, with the strength constrained by 
the trap spacing (J <C h(ü osc ). The relative strength of the in- 
teractions along z and x can be changed by rotating the angle 
Ç of polarization out of plane. Interactions between more dis- 
tant neighbors are relatively weak because the far off resonant 
coupling at larger r cannot distinguish the spin dependence of 
excited states. 




1 2 3 4 2 4 
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FIG. 3: (Color online) Design and verification of noise protected 
ground states arising from a simulation of H^ ín . The system is com- 
prised of 9 polar molècules trapped in a 3 x 3-square-lattice in the 
z — f plane with lattice spacing b = ry/ ' \fï driven with a field of fre- 
quency (B/r and out of plane polarization angle Ç. (a) Noise resilience 
of the ground states of the resultant interaction H sp [ a as a function 
of field frequency, quantified by the rms-magnetisations of the two 
ground-states, 8 2 S Z = 8 2 S A (solid-lines) and 8 2 S- V (dashed lines) for 
Ç = 0. The inset shows the protected region, when tuning near to 
the 2g resonance E(2 g ) ~ 1.9y which realizes the model üf^. (b) 
Absorption-spectroscopy of ground states /"(c^probe) for two spin- 
textures obtained by tuning C0/r near the 2 g resonance, 8/y= 1.88, 
with b = ry/ \/2. For Ç = the spectrum is gapped by J /2, which 
is a signature of a protected qubit (top), while for Ç = ir/2 the ex- 
citations are gap-less spin-waves (bottom). The peak at co plo b e = 
is of order unity. The spectroscopic coupling component a = x,y,z 
(solid,dashed,dotted lines) are indicated. 



The authors of Ref. [8] show that the ideal spin model 
I (for Ç 7^ ±7r/2) has a 2-fold degenerate ground subspace, 
which is gapped with weak size dependence for cos Ç = 1 . 
The the two ground-states, which we denote, |0)l and 
have zero local magnetizations (of } )L- Our implementation is 
not ideal because there are residual longer range interactions 
along several directions in the plane as well as off resonant 
couplings to excited state potentials yielding unwanted spin 
patterns. We note, however, that all the effective spin Hamil- 
tonians described in Eq. (5) obtained using fields with linear 
polarization involve sums of products of Pauli operators and 
hence are invariant under time-reversal. For l odd, the degen- 
eracy present in the ground state of H^ a is of Kramers' type 
and imperfect implementation will not break this degeneracy 
though it may decrease the energy gap. 

We have numerically computed the effective interaction on 
a £ 2 — 3 x 3 square lattice with spacings b = ry/ V2 and we 
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take the localization to the point dipole limit. In Fig. 3(a) we 
plot the oc = x,y,z-components of the rms magnetization for 
the ground subspace, 8 2 S a = ZijZao \l(G'\o?j\G) l \ 2 /2f, 
as a function of the detuning cop — 2Bj li for polarization an- 
gle Ç = 0. This allows for computing logical qubit errors 
due to quasi-static noise. Near the bare resonance h(üf — 
2B = y/2 the system show múltiple long-range resonances 
as all the sites couple near-resonantly at coupling strength 
°c l/£> 3 . The last of these long-range resonance appears at 
/icOf -2Bk 1.36y for the interaction between next nearest 
neighbor sites with spacings of \fï.b. The 2 g -resonance lies 
at h(úf —2Brí 1.9y for nearest neighbor sites and shows the 
remarkable feature of no magnetization on any site in any 
space-direction oc withing the ground-state manifold (see in- 
set). The resulting immunity of the system to local noise can 
be probed by applying an homogeneous fi-field of frequency 
C0p ro b e polarized in the direction oc = x,y,z- The correspond- 
ing absorption spectrum for an arbitrary code state \\\t)i is, 
X a (co prob e) = -rir3{ L (y\S a (r«a pmbe -H svin + ihr)- l S a \y) L } 
where S a = LijOfj/f and F is an effective linewidth. This 
quantity is plotted in Fig. 3(b) for two different spin-textures 
obtained for the same field frequency C0/r = 1 .88y but different 
polarizations and were we set F = 0. 1/. For polarization Ç = 
(see top inset) one realizes the protected qubit, whose spec- 
trum is gapped by J/2. For polarization along the i-direction 
Ç = jt/2 (see bottom inset) the ground-subspace is given by a 
set of l quantum-Ising stripes along z, whose spectrum is un- 
gapped with a large peak at CO pr obe = in response to coupling 
with a B field polarized along oc = x. 



Spin model II is likewise obtained using this mechanism. 
Consider a system of four molècules connected by three 
length b edges forming an orthogonal triad in space. There 
are several different microwave field configurations that can 
be used to realize the interaction H^ n along the links. One 
choice is to use two microwave fields polarized along z, 
one tuned near resonance with a 1 ? potential and one near 
a 1, potential. When the detunings and Rabi freqüències 
are chosen so that (|OiJC(l g ,3,3) - |QijC(l,„ 1, l)) re i = 
then the resultant spin pattern is Eq. (2) with Jj_ = 
-ft<|ni f |C(l„3,3)W4 and = »|<|n lí |C(l í ,2,2)) B a/4. 
The ratio |/i|/|/ z | can be tuned either by changing the lat- 
tice spacing or by using a third microwave field polarized 
along z and tuned near the 2 g potential, in which case J± — » 
J± +h(\Q.2 g \C(2 g )} Is i/8, A bipartite lattice composed of such 
triads with equally spaced nearest neighbors can be built using 
two planes of stacked triangular lattices. Such a lattice could 
be designed using bichromatic trapping làsers in two spatial 
dimensions and a suitably modulated lattice in the third di- 
mension normal to both planes. A realization of model II us- 
ing a different set of 3 microwave fields is shown in Fig. 4. 
The obtained interaction is close to ideal with small residual 
coupling to next nearest neighbors as in model I. 




FIG. 4: (Color online) Implementation of spin model ^, Shown 
is the spatial configuration of 12 polar molècules trapped by two 
parallel triangular lattices (indicated by shaded planes) with sepa- 
ration normal to the plane of Az/v3 and in plane relative lattice shift 
of Azy/2/3. Nearest neighbors are separated by b = Az and next 
nearest neighbor couplings are at \Í2b. The graph vertices repre- 
sent spins and the edges correspond to pairwise spin couplings. The 
edge color indicates the nature of the dominant pairwise coupling for 
that edge (blue= a z a z , red= a y a y , green= a v a l , black="other"). 
For nearest neighbor couplings, the edge width indicates the relative 
strength of the absolute value of the coupling. For this implemen- 
tation, the nearest neighbor separation is b = ry. Three fields all 
polarized along z were used to generate the effective spin-spin in- 
teraction with freqüències and intensities optimized to approximate 
the ideal model H^ a . The field detunings at the nearest neighbor 

spacing are: h(úi -£(l g (l/2)) = -0.05y/2, to 2 - E(0 g (1/2)) = 
0.05y/2,to 3 -£(2g(3/2)) =0.10y/2and the amplitudes are |£2i| = 
4|ÍÍ2| = |ÍÏ3| = O.Oly/A. For y= 40MHz this generates effective 
coupling strengths J z = — 100kHz and = —0.4./-. The magni- 
tude of residual nearest neighbor couplings are less than 0.04|/ ; | 
along x and y-links and less than 0.003|/-| along z-links. The size of 
longer range couplings J\ r are indicated by edge line style (dashed: 
|/lr] < 0.01|7,|, dotted: |y lr | < 10" 3 |y ; |). Treating pairs of spins 
on z-links as a single effective spin in the low energy sector, the 
model approximates Kitaev's 4-local Hamiltonian [10] on a square 
grid (shown here are one palquette on the square lattice and a neigh- 
bor plaquette on the dual lattice) with an effective coupling strength 
4ff = ~(J± M:)Vzl/16 ~ 167Hz. 



VI. CONCLUSIONS 

We have shown how to engineer pairwise spin- 1/2 (or qubit) 
interactions between polar molècules with designable range 
and spatial anisotropy. The couplings are strong relative to 
decoherence rates, and the simulation does not require com- 
plicated control sequences to build effective time averaged 
Hamiltonians, spin-dependent lattices or local addressability. 
Thus polar molècules in optical lattices are a good candidate 
to provide a laboratory realization of exòtic states of matter. 
We envision that engineering of these new materials might 
eventually provide the basis of a new type of quantum òp- 
tics, where systematic tools are developed which allow the 
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controlled preparation and manipulation of excitations such 
as anyons, with applications ranging from fundamental qües- 
tions in condensed matter physics to quantum computing. 
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APPENDIX A: EFFECTIVE INTERACTIONS 



Here the energy dependent terms í(|y|J(/)) = /z|£2|/[7ïCOf — 
£■(171^(7))] quantify the amplitude in the excited states. 
The energies E(\Y\^(J))) correspond to eigenvalues of H mt 
and the sets {Kj(Q%(J))} 2 J=1 and {KJ( a l a (J))} 3 j=l are co- 
efficients of the eigenvectors for |F| = 0,1. For Y = 
0, the energies are £(0+(l ±1/2)) = 2B + y[a3d 2 /2yr 3 - 
1 /4 ± y/(ad 2 /2yr 3 + 1 /4) 2 + 1/2] , and £(0 o (l ± 1/2)) = 
2B+j[-ad 2 /2yr 3 - 



•l/4) 2 + 1/2 _ 
1/4 + y/(-c3d 2 /2yr 3 - 



1/4)2+1/2], 

The eigenvector components are K\(0q) = cos(Yo»>/2) and 
•^2(0™) = sin(y «i/2) where the angles satisfy tan(y +) = 
v%/ '(1/2 + ad 2 /yr 3 ), and tan(y 0o ) = V2/{l/2-a3d 2 /yr 3 ). 
For Y = ±1 the eigenvectors and doubly degenerate eigenval- 
ues are obtained by diagonalizing the 3x3 matrices: 



The effective spin-spin interaction Eq. (5) between polar 
molècules depends both on the frequency C0/r and polarization 
e/r = oc_ è- 1 + ao<?o + oc+é 1 , (éo = 2) of the field. The explicit 
form for the coupling coefficients is: 

A u = |ao| 2 [C(0-,l,2)-C(0+ 1,2)] 

+(|a_| 2 + |a + | 2 )[C(l i ,,3,3)-C(l 1 ,,l,l)] 
+9t[<a_][C(0-,2,l)-C(0+ 2,1)], 

A 2 ,2 = A M -29l[o^a-][C(Oj,2,l)-C(0+2,l)], 

A 3 , 3 = |a | 2 [2C(l„2,2)-C(0 g ,l,2)-C(0+ 1,2)] 
+(|a + | 2 + |a_| 2 )[C(2,) + C(0+,2,l)/2 
+C(0j, 2,l)/2-C(l Ul l,l)-C(l gl 3,3)], 

Ai, 2 = 3[a+a_](C(0-,2,l)-C(0+ 2,1)), 

Ai j3 = 9t[a + a -a5a_]C(l 4 ,2,3), 

A 2 , 3 = 3[<oco-oc5a-]C(l g ,2,3), 

A ,i = 9ï[a+a + a5a_]C(l i? ,2,3), 

Aq.2 = 3[a^oco + oc5a-]C(l g ,2,3), 



A ,3 



a 4 



|a_| 2 )[C(2,)-C(0+ 2,l)/2 



-C(0",2,l)/2]. 



The component Ao,o weights a scalar energy shift which we 
ignore. The coefficients C(|F|J ) quantify coupling to excited 
states with different symmetries and are given by 

C(0»J,*) = Kj(0'») 2 s(0<S(3/2))+K k (0>S) 2 s(0>»(l/2)), 
C(lc,j,k) = lti^("la(3/2))^( fl l a (3/2)),(«l (3/2)) 

+ LLl^( Ò la(l/2))^( è la(l/2)),( è lo(l/2)) 

C(2 g ) = s{2 g (3/2)). 



2B1 3 



í-a 2 4 

±1 



±1 

-a 4 4 



1 

yr 1 / 



For Y = ±2, the eigenvalues are doubly degenerate with ener- 
gies £(2 (3/2)) = 2B + y/2 + ad 2 /r 3 . 

A caveat is that we do not have point dipoles but rather 
wavepackets with spatial distributions parallel and perpen- 
dicular to the intermolecular axis f. Components of inter- 
molecular separations orthogonal to z will couple to states 
with different symmetry and an exact treatment would re- 
quire averaging over the angular distrubution with the ap- 
propriate frame transformation. However, we argue that in 
our regime this finite size affect is negligible. The relative 
magnitude can be estimated by the ratio of the marginal rel- 
ative coordinate probability distributions perpendicular and 
parallel to £. Defining p±(r) = J í/í2sin 2 0r 2 |\|/ re i(r,0)| 2 and 
p\\(r) = J c/í2cos 2 0r 2 |\)/ r ei(r,9)| 2 , the peak of the distribu- 
tions is at r = Az where for zo/Az *C 1, the relative amount 
of unwanted couplings is pj_(Az)/ p^(Az) ~ 4(zo/Az) 2 . For 
molecular wavepacket localization 27tzoAtrap = 0.1, the ratio 
is pj_(Atop)/p||(Àtrap) ~ 10 -3 , hence it is warranted to com- 
pute the couplings as if the entire weight of the wavefunction 
were parallel to £. 
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